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Section  I 


Introduction 

In  Reference  1,  we  develop  the  probability  distribution  for  demand  in  a  lead  time 
for  the  situation  in  which  customers  arrive  according  to  a  Poisson  process  and  each 
customer  requires  X  units,  where  X  is  described  by  the  logarithmic  distribution;  i.e.  the 
probability  mass  function  of  X  is 

f  (x)  =  "x  iKe"  x  =  l’  2»  ‘  *  * 

In  addition,  lead  times  are  assumed  to  be  Gamma  Distributed.  We  call  the  resulting 
distribution  of  demand  in  a  lead  time  the  Logarithmic-Poisson-Gamma  (LPG)  distribu¬ 
tion.  This  distribution  appears  to  be  a  useful  approximation  to  the  demand  processes  of 
a  number  of  Air  Force  EOQ-type  (consumable)  items. 

Analytic  details  of  the  LPG  distribution  are  presented  in  Reference  I.  In  addition, 
reference  1  presents  a  scaled  negative  binomial  approximation  to  the  LPG.  We  assume 
the  reader  is  familiar  with  these  results.  In  this  paper,  we  present  a  sensitivity  analysis 
of  the  LPG  distribution  and  of  the  associated  scaled  negative  binomial  approximation. 


SECTION  II 


Sensitivity  Analysis 

To  get  a  feel  for  the  shape  of  the  LPG  distribution,  and  to  evaluate  the  accuracy  of 
the  scaled  negative  Binomial  approximation,  we  computed  the  cumulative  distribution 
function  values  for  both  of  these  distributions  for  40  different  parameter  sets.  The 
BASIC  computer  code  used  for  these  calculations  is  presented  in  Appendix  B,  while 
Figure  II- 1  illustrates  the  results  of  our  calculations  for  one  of  these  data  sets.  This 
figure  presents  the  cumulative  distribution  function  (CDF)  and  probability  density 
function  (PDF)  for  both  the  LPG  and  scaled  negative  binomial  distributions  for  the  case 
in  which  the  average  requisition  size  is  2  units,  the  average  monthly  demand  rate  is  .5, 
the  average  lead  time  is  8  months,  and  the  coefficient  of  variation  of  lead  time  is  .25. 
For  this  situation,  the  mean,  coefficient  of  variation,  skewness,  and  kurtosis  of  the 
distribution  of  demand  in  a  leadtime  equals  4.00,  0.97,  1.46,  and  6.40,  respectively. 
These  values  are  shown  on  the  left  hand  side  of  Figure  II- 1.  The  specific  parameters 
utilized  in  the  LPG  distribution  calculations  are  also  shown  on  the  left  hand  side  of  this 
figure,  as  well  as  specific  CDF  values.  The  column  labeled  "FX"  presents  cumulative 
distribution  function  values  for  the  LPG  distribution,  while  the  column  labeled  "FNBX" 
presents  the  corresponding  cumulative  distribution  function  values  for  the  scaled 
negative  binomial  approximation.  Finally,  the  column  labeled  "DIFF"  presents  the 
difference  between  the  LPG  distribution  values  and  the  corresponding  negative  binomial 
CDF  values. 

The  right  hand  side  of  Figure  II- 1  plots  both  the  cumulative  and  density  function 
values  associated  with  the  LPG  and  scaled  negative  binomial  distributions.  As  shown  in 


» 


n-3 


the  Figure,  the  scaled  negative  binomial  approximation  underestimates  the  correspond¬ 
ing  LPG  CDF  values  for  low  levels  of  demand  in  the  leadtime,  but  the  approximation 
becomes  much  better  for  demand  values  exceeding  10  units.  Note  the  slight  down-turn 
labeled  as  Point  A.  This  is  caused  by  the  rounding  rules  employed  in  the  scaled  negative 
binomial  code.  A  slight  change  in  these  rules  may  be  used  to  make  the  CDF  monotonic. 

Figure  II-2  presents  a  similar  set  of  calculations  in  which  the  monthly  demand  rate 
is  equal  to  2  units  per  month,  and  the  other  parameters  are  the  same  as  used  in  Figure 
11-1.  In  this  case,  the  scaled  negative  binomial  distribution  is  a  much  closer 
approximation  to  the  LPG  curve. 

We  will  refer  to  the  specific  data  sets  plotted  in  Figures  II- 1  and  II-2  as  the  "LO" 
and  "HI"  base  cases,  respectively.  Table  II- 1  presents  the  specific  parameter  values 
selected  for  these  calculations.  The  data  set  numbers  assigned  to  each  of  these  curves 
are  shown  in  the  right  hand  side  of  Table  II- 1  and  detailed  plots  associated  with  each  of 
the  data  sets  are  presented  in  Appendix  A.  Let  us  now  discuss  our  results. 

Results  for  LO  Base  Case. 

As  noted  above,  the  LO  Base  Case  has  the  following  characteristics: 

Average  Requisition  Size  (Units) 

Average  Units  Demanded  per  Month 
Average  Lead  Time  (Months) 

Coefficient  of  Variation  of  Lead  Time 


Value 
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.5 
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Figure  II-2.  LPG  and  Scaled  Negative  Binomial  Distributions  for  the  HI  Demand  Base  .Case 


1 1-5 


Demand  Rate  in  .1, .5, 1., 1.5, 2.0 

Units  per  Month 

Coefficient  of  Variation  .01,  .25,  .50,  .75,  1.0 
of  Leadtime 
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Let  us  now  consider  the  impact  of  variations  in  these  characteristics. 

Requisition  Size  Sensitivity.  In  performing  our  sensitivity  analyses,  three  of  the 
major  distribution  characteristics  defined  in  Table  11-1  were  held  constant,  while  the 
fourth  characteristic  was  varied.  For  example,  data  set  numbers  1  thru  5  were 
constructed  by  holding  average  leadtime,  demand  rate,  and  the  coefficient  of  variation 
of  leadtime  equal  to  the  values  for  the  LO  demand  base  case,  while  varying  the  average 
requisition  size.  Data  set  numbers  21  thru  25  were  performed  similarly,  but  in  this  case 
the  HI  demand  base  case  was  used. 

Details  of  our  results  fo t  all  40  data  sets  are  presented  in  Appendix  A,  while  Table 
II-2  and  Figures  U-3  thru  II-7  illustrate  our  results.  Table  II-2  summarizes  the  primary 
calculations  associated  with  data  sets  1  thru  5,  while  Figures  II-3  thru  II-7  present  the 
associated  plots  of  the  LPG  CDF.  The  table  summarizes  the  sensitivity  of  the  primary 
LPG  distribution  parameters  (i.e.  Theta  Lambda,  Alpha,  and  Beta)  to  changes  in  average 
requisition  size  for  the  low  demand  case  base.  This  table  displays  the  requisition  size, 
demand  rate,  mean  leadtime  and  leadtime  variability  coefficients  that  were  used  to 
drive  the  data  set  calculations,  as  well  as  the  associated  LPG  parameters  and  the 
moments  of  the  LPG  leadtime  demand  distribution.  Finally,  at  the  bottom  of  Table  II- 
2,  we  present  the  specific  demand  values  associated  with  given  percentage  points  of  the 
LPG  distribution  and  of  the  associated  scaled  negative  binomial  approximation.  For 
example,  consider  the  data  shown  in  the  first  column  at  the  bottom  of  Table  U-2.  For 
data  set  number  1,  a  reorder  point  of  eight  units  is  required  to  insure  a  95%  probability 
of  no  stockouts  during  a  leadtime.  Using  the  scaled  negative  binomial  approximation 
for  data  set  1,  a  reorder  level  of  9  units  is  required  to  produce  the  same  93% 
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SAMPLE  LPG  CALCULATIONS 


Data  Set  No. 

1 

2 

3 

4 

5 

Ave.  Req.  Size 

1.01 

2 

4 

8 

16 

Ave.  Demand /Mo. 

.5 

.5 

.5 

.5 

.5 

Mean  Lead  Time 

8.0 

8.0 

8.0 

8 

8 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 
Theta 

.02 

.72 

.90 

.96 

.99 

Lambda 

.49 

.25 

.12 

.06 

.03 

Alpha 

16 

16 

16 

16 

16 

Beta 

2 

2 

2 

2 

2 

LPG  Moments 

Mean 

4 

4 

4 

4 

4 

CV 

.56 

.97 

1.63 

2 

.63 

4.15 

Sk ewness 

.41 

1.46 

2.96 

5 

.08 

8.18 

Kurtosis 

3.23 

6.40 

16.45 

42 

.13 

103.79 

1 

2 

3 

4 

5 

Percentage  LPG 

Points 

NB 

fPCf 

El 

LPG 

El 

lPg 

""Eb 

Lpg 

NB 

.50  4 

5 

3 

4 

1 

4 

0 

8 

0 

18 

.60  4 

5 

4 

5 

2 

5 

0 

10 

0 

22 

.70  5 

.  6 

5 

6 

4 

6 

2 

11 

0 

26 

.80  6 

7 

7 

7 

7 

10 

4 

13 

1 

29 

.85  6 

7 

8 

8 

9 

12 

7 

14 

3 

31 

.90  7 

8 

9 

9 

12 

14 

12 

22 

7 

33 

.95  8 

9 

12 

12 

17 

16 

22 

33 

22 

43 

.97  9 

10 

13 

13 

21 

22 

31 

38 

38 

68 

.99  10 

11 

17 

16 

31 

26 

52 

42 

82 

93 

Table  II  2.  Requisition 

Size  Sensitivity 

for  LO 

Base 

Case 
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probability.  Scanning  the  results  for  all  five  data  sets,  we  observe  that  the  scaled 
negative  binomial  provides  reorder  point  results  very  similar  to  the  LPG  distribution  for 
average  requisition  sizes  of  4  or  less,  but  provides  a  very  poor  approximation  for  the 
high  average  requisition  size  values. 

Lead  Time  Sensitivity.  Table  II-3  presents  percentage  points  for  the  LPG  and 
Scaled  Negative  Binomial  (SNB)  approximations  for  the  LO  Base  case  and  Data  Sets  6 
thru  10.  This  Table  shows  the  effect  of  mean  lead  times  of  4,  6,  S,  10,  and  12  months, 
respectively.  Note  that  the  SBN  yields  percentage  points  which  are  very  close  to  the 
LPG  values,  and  that  the  approximation  improves  as  the  lead  time  increases. 

Demand  Rate  Sensitivity.  Table  II-4  presents  results  of  varying  the  LO  Base  case 
using  demand  rates  of  .1,  .5,  1.0,  1.5,  and  2.0  units  per  month,  respectively,  Note  that 
the  LPG  and  SNB  values  are  vary  similar  for  all  but  the  demand  =.l  case. 

Lead  Time  Variability  Sensitivity.  Table  II-5  presents  the  results  of  increasing  lead 
time  variability  (Dat  Sets  16-20).  Observe  the  percentage  points  are  within  1  unit  for 
all  but  the  C V=.75,  60%  point.  Hence,  the  approximation  is  excellent  for  all  five  data 
sets. 


Results  for  the  HI  Base  Case.  The  HI  Base  Case  has  the  following  characteristics: 

Value 

Average  Requisition  Size  2 

Average  Units  Demanded  per  Month  2 

Average  Lead  Time  (Months)  8 

Coeeficient  of  Variation  of  Leat  Time 


.25 
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Data  Set  No. 

SAMPLE 

6 

LPG  CALCULATIONS 

7  8 

_ 9 

10 

Ave.  Req.  Size 

2 

.00 

2.00 

2 

.00 

2.00 

2 

.00 

Ave.  Demand /Mo 

.50 

.50 

.50 

.50 

.50 

Mean  Lead  Time 

4, 

.00 

6.00 

8 

.00 

10.00 

12 

.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.25 

.25 

.25 

.25 

.25 

Alpha 

16 

.00 

16.00 

16 

.00 

16.00 

16 

.00 

Beta 

4, 

.00 

2.67 

2 

.00 

1.60 

1 

.33 

LPG  Moments 

Mean 

2, 

.00 

3.00 

4 

.00 

5.00 

6 

.00 

Skewness 

2, 

.16 

1.72 

1 

.46 

1.28 

1 

.14 

Kurtosis 

10 

.22 

7.67 

6 

.40 

5.65 

5 

.15 

6 

7 

8 

9 

10 

Percentage  LPG 

HE 

LPG- 

nB 

LPG 

El 

Lt>G 

HE 

LPG 

NB 

Points 

.50  1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

.60  2 

2 

3 

4 

4 

5 

5 

6 

6 

7 

.70  2 

4 

4 

5 

5 

6 

6 

7 

8 

8 

.80  3 

5 

5 

6 

7 

7 

8 

9 

9 

11 

.85  4 

5 

6 

7 

8 

8 

9 

9 

11 

12 

.90  5 

6 

7 

8 

9 

9 

11 

12 

12 

13 

.95  7 

8 

10 

9 

12 

12 

14 

14 

15 

15 

.97  9 

9 

11 

12 

13 

13 

15 

15 

17 

16 

.99  12 

9 

15 

14 

17 

16 

19 

19 

22 

20 

Table  1 1 -3.  Lead  Tine  Sensitivity  for  LO  Base  Case 
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SAMPLE 

LPG  CALCULATIONS 

Data  Set  No. 

11 

12 

13 

14 

15 

Ave.  Req.  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand/Mo 

.10 

.50 

1.00 

1.50 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.05 

.25 

.50 

.75 

1.00 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2.00 

2.00 

2.00 

LPG  Moments 

Mean 

.80 

4.00 

8.00 

12.00 

16.00 

CV 

2.11 

.97 

.71 

.60 

.53 

Skewness 

3.52 

1.46 

.95 

.73 

.60 

Kurtosis 

21.73 

6.40 

4.53 

3.93 

3.64 

11 

12 

13 

14  15 

Percentage  LPG 

Points 

NB  CPC 

NB 

Tpg  Ef 

LPG 

~ CB  LPG 

""nB 

.50 

0 

1  3 

4 

7  8 

11 

12  15 

15 

.60 

0 

2  4 

5 

8  9 

13 

13  17 

18 

.70 

1 

2  5 

6 

10  11 

15 

15  19 

20 

.80 

1 

2  7 

7 

12  13 

17 

19  23 

23 

.85 

2 

2  8 

8 

14  14 

19 

20  26 

26 

.90 

3 

5  9 

9 

16  16 

22 

22  27 

28 

.95 

4 

5  12 

12 

19  19 

25 

26  32 

32 

.97 

5 

6  13 

13 

21  21 

28 

28  35 

35 

.99 

8 

8  17 

16 

26  25 

34 

33  41 

40 

Table  II *4.  Demand  Rate  Sensitivity  for  LO  Base  Case 
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SAMPLE 

LPG  CALCULATION 

Data  Set  No. 

16 

17 

18 

19 

20 

Ave.  Req. 

Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand /Mo. 

0.5 

0.5 

0.5 

0.5 

0.5 

Mean  Lead 

Time 

8 

8.0 

8.0 

8.0 

8.0 

CV  of  Lead 

Time 

.01 

.25 

.50 

.75 

1.0 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

• 

.25 

.25 

.25 

.25 

.25 

Alpha 

1 

,000 

16.00 

4.00 

1.78 

1.00 

Beta 

1 

,250 

2.00 

.50 

.22 

.13 

LPG  Moments 

Mean 

4.00 

4.00 

4.00 

4.00 

4.00 

CV 

.94 

.97 

1.06 

1.20 

1.37 

Skewness 

1.61 

1.46 

1.21 

1.13 

1.29 

Kurtosis 

6.84 

6.40 

5.55 

5.15 

5.76 

Percentage 

Points 

16 

17 

18 

19  20 

IFg~ 

“TO 

LPG 

“TO 

LEG  nb 

LPG 

TO  IE! 

"TO 

.50 

3 

4 

3 

4 

3  4 

2 

3  2 

2 

.60 

4 

5 

4 

5 

4  5 

3 

5  3 

4 

.70 

5 

6 

5 

6 

5  6 

5 

6  5 

6 

.80 

7 

7 

7 

7 

7  8 

7 

8  7 

8 

.85 

8 

8 

8 

8 

8  9 

8 

9  9 

9 

.90 

9 

9 

9 

9 

10  11 

10 

11  11 

12 

.95 

11 

12 

12 

12 

12  13 

14 

14  15 

15 

.97 

13 

13 

13 

13 

14  14 

16 

17  18 

19 

.99 

17 

15 

17 

16 

19  19 

21 

21  25 

25 

Table  1 1 -5.  Lead  Tine  Variability  Sensitivity  for  LO  Base  Case 


i 


Recall  that  Figure  II-2  plots  the  LPG  and  SNB  curves  for  this  case,  and  that  the 
detailed  plots  for  each  data  set  are  presented  in  Appendix  A. 


Requisition  Size  Sensitivity.  Table  11-6  presents  the  results  for  Data  Sets,  21-25, 
reflecting  average  requisition  sizes  of  1.01,  2.0,  4.0,  8.0,  and  16.0,  respectively.  Note 
that  the  LPG  and  SNB  curves  yield  similar  results  for  average  requisition  sizes  of  4  or 
less,  but  differ  significantly  for  average  requisition  sizes  of  8  and  16. 

Lead  Time  Sensitivity.  Table  II-7  presents  results  for  Data  Sets  26-30,  reflecting 
lead  times  of  4,  6,  8,  10,  and  12  months,  respectively.  For  this  case,  the  LPG  and  SNB 
yield  similar  results  in  ail  cases. 

Demand  Rate  Sensitivity.  The  results  for  Data  Sets  31-35  are  presented  in  Table 
II-8.  These  Data  Sets  have  average  demands  of  .1,  .5,  1,  1.5,  and  2  units  per  month, 
respectively.  Observe  that  the  SNB  approximation  is  a  particularly  good  fit  to  the  LPG 
except  for  the  50  and  60%  points  for  the  .1  units  per  month  case. 

Lead  Time  Variability  Sensitivity.  Table  II-9  presents  results  for  coefficients  of 
variation  of  .01,  .25,  .50,  .75,  and  1.0,  respectively  (Data  Sets  36-40).  Note  that  the 
LPG  and  SNB  distributions  are  very  similar  for  all  of  these  cases. 


SAMPLE  LPG  CALCULATIONS 


Data  Set  No. 

21 

22 

23 

24 

25 

Ave.  Req.  Size 

1.01 

2.00 

4.00 

8.00 

16.00 

Ave.  Demand /Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.02 

.72 

.90 

.96 

.99 

Lambda 

1.98 

1.00 

.50 

.25 

.13 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2.00 

2.00 

2.00 

LPG  Moments 

Mean 

16.00 

16.00 

16.00 

16.00 

16.00 

CV 

.36 

.53 

.84 

1.33 

2.09 

Skewness 

.27 

.60 

1.35 

2.44 

4.02 

Kurtosis 

3.13 

3.64 

5.97 

12.32 

27.79 

21 

22 

23 

24  25 

Percentage 

Joints 

IFg~ 

NB  L^G 

LPG  NB 

IPg 

"175  LPd" 

.50 

15 

16  15 

15 

13  14 

8 

13  2 

22 

.60 

17 

18  17 

18 

16  19 

12 

18  5 

26 

.70 

19 

19  19 

20 

20  23 

18 

26  11 

31 

.80 

21 

21  23 

23 

25  27 

27 

35  23 

37 

.85 

22 

22  25 

26 

29  31 

33 

39  33 

56 

.90 

24 

24  27 

28 

34  35 

43 

47  50 

75 

.95 

26 

27  32 

32 

42  42 

59 

61  80 

93 

.97 

28 

28  35 

34 

49  46 

72 

67  105 

101 

.99 

31 

32  41 

40 

61  57 

99 

88  163 

148 

Table  II -6.  Requisition  Size  Sensitivity  for  HI  Base  Case 
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SAMPLE 

LPG  CALCULATIONS 

Data  Set  No. 

26 

27 

28 

29 

30 

Ave.  Req  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand /Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead  Time 

4.00 

6.00 

8.00 

10.00 

12.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

1.00 

1.00 

1.00 

1.00 

1.00 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

4.00 

2.67 

2.00 

1.60 

1.33 

LPG  Moments 

Mean 

8.00 

12.00 

16.00 

20.00 

24.00 

CV 

.71 

.60 

.53 

.49 

.46 

Skewness 

.95 

.73 

.60 

.52 

.47 

Kurtosis 

4.53 

3.92 

3.64 

3.48 

3.38 

26 

27 

28 

29  30 

Percentage 

Points 

LPG 

“Tif  TFT 

“KB 

LPG  NB 

LPG 

TJB  LPG 

~R| 

.50 

7 

8  11 

12 

15  15 

19 

20  23 

23 

.60 

8 

9  13 

13 

17  18 

21 

22  25 

26 

.70 

10 

11  15 

15 

19  20 

24 

25  29 

29 

.80 

12 

13  17 

19 

23  23 

28 

28  33 

33 

.85 

14 

14  19 

20 

25  26 

30 

30  35 

36 

.90 

16 

16  22 

22 

27  28 

33 

34  39 

39 

.95 

19 

19  25 

26 

32  32 

38 

37  44 

44 

.97 

21 

21  28 

28 

35  34 

41 

41  48 

48 

.99 

26 

25  34 

33 

41  40 

48 

48  55 

55 

Table 

II  "7 

.  Lead  Time  Sensitivity  for 

HI  Base  Case 
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SAMPLE 

LPG  CALCULATIONS 

Data  Set  No, 

31 

32 

33 

34 

35 

Ave.  Req  Size 

2.00 

2.00 

2 

.00 

2.00 

2.00 

Ave.  Demand /Mo. 

.10 

.50 

1 

.00 

1.50 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8 

.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.05 

.25 

.50 

.75 

1.00 

Alpha 

16.00 

16.00 

16 

.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2 

.00 

2.00 

2.00 

LPG  Moments 

Mean 

.80 

4.00 

8 

.00 

12.00 

16.00 

CV 

2.11 

.97 

.71 

.60 

.53 

Skewness 

3.52 

1.46 

.95 

.73 

.60 

Kurtosos 

21.73 

6.40 

4 

.53 

3.92 

3.64 

31 

32 

33 

34  35 

Percentage  LPG 

Points 

Nb  CPG 

“"TF 

Lpg 

“ TJF 

LP(S 

IE? 

.50 

0 

1  3 

4 

7 

8 

11 

12  15 

15 

.60 

0 

2  4 

5 

8 

9 

13 

13  17 

18 

.70 

1 

2  5 

6 

10 

11 

15 

15  19 

20 

.80 

1 

2  7 

7 

12 

13 

17 

19  23 

23 

.85 

2 

2  8 

8 

14 

14 

19 

20  25 

26 

.90 

3 

5  1 

9  9 

16 

16 

22 

22  27 

28 

.95 

4 

5  12 

12 

19 

19 

25 

26  32 

32 

.97 

5 

6  13 

3 

21 

21 

28 

28  35 

34 

.99 

8 

8  17 

16 

26 

25 

34 

33  41 

40 

Table  11*8.  Demand  Rate  Sensitivity  for  HI  Base  Case 


11-20 


SAMPLE 

LPG  CALCULATIONS 

Date  Set  No. 

36 

37 

38 

39 

40 

Ave.  Req.  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Pemand/Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CY  of  Lead  Time 

.01 

.25 

.50 

.75 

1.00 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

1.00 

1.00 

1.00 

1.00 

1.00 

Alpha 

10,000 

16.00 

4.00 

1 . 78 

1.00 

Beta 

1,250 

2.00 

.50 

.22 

.13 

LPG  Moments 

Mean 

16.00 

16.00 

16.00 

16.00 

16.00 

CV 

.47 

.53 

.69 

.88 

1.10 

SVewness 

.80 

.60 

.65 

1.03 

1.55 

Kurtosis 

3.96 

3.64 

3.69 

4.86 

7.10 

36 

37 

38 

39  40 

Percentage 

LPG 

NB 

U>G 

KB 

Lpg  Kb 

rpG 

™KB  LPG 

“KB 

Points 

.50 

15 

16 

15 

15 

14  14 

12 

13  10 

11 

.60 

17 

18 

17 

18 

17  18 

16 

16  14 

15 

.70 

19 

20 

19 

20 

20  20 

19 

20 

.80 

22 

22 

23 

23 

24  25 

26 

23  27 

27 

.85 

24 

25 

25 

26 

27  27 

30 

26  32 

33 

.90 

26 

27 

27 

28 

31  31 

35 

35  39 

40 

.95 

30 

29 

32 

32 

37  37 

44 

44  51 

51 

.97 

32 

32 

35 

34 

41  41 

50 

50  60 

61 

.99 

38 

36 

41 

40 

50  51 

64 

63  80 

80 

Table  II -9. 

Lead 

Time  Variability  Sensitivity  for 

HI  Base  Case 

II 


Moments  of  the  LPG  Distribution 


As  noted  above,  we  computed  the  first  four  moments  of  the  LPG  distribution. 
Figures  11-8  and  II-9  present  the  means  and  standard  deviations  associated  with  each  of 
the  forty  data  sets  studied,  while  Figure  H-10  presents  a  plot  of  (skewness)^  value 
versus  the  kurtosis  for  the  forty  data  set  calculations.  Figure  H-10  is  particularly 
interesting.  As  shown  in  Johnson  and  Kotz  (1969,  Pg  39),  the  skewness^-kurtosis 
relationships  of  many  well  known  probability  distributions  are  described  by  specific 
lines,  points,  or  regions  when  plotted  on  a  chart  similar  to  Figure  H-10.  For  example, 
the  lines  and  regions  associated  with  Normal,  Gamma,  Log-normal,  Poisson,  Binomial, 
and  Negative  Binomial  Distributions  are  shown  in  Figure  II- 10.  Hence,  the  plot  of 
skewness^  versus  kurtosis  for  a  specific  probability  distribution  may  be  used  to  gain 
insights  into  appropriate  approximations  for  that  distribution.  As  shown  in  Figures  H-10 
thru  11-13,  all  of  the  sample  skewness^-kurtosis  points  associated  with  our  LPG 
calculations  fall  on  a  line  which  is  almost  identical  with  that  associated  with  the 
gamma  probability  distribution.  Hence,  it  appears  that  the  gamma  may  be  a  very  useful 
approximation  for  LPG  calculations. 


Standard  Deviation. 


Section  III 


Summary  and  Conclusions 

In  this  paper,  we  have  presented  sample  distribution  functions  for  the  Logarithmic- 
Poisson-Gamma  (LPG)  distributed  and  for  the  Scaled  Negative  Binomial  (SNB)  approxi¬ 
mation  to  the  LPG.  Two  base  cases  were  considered;  these  cases  have  the  following 
characteristics. 


LO 

HI 

Base 

Base 

Case 

Case 

Average  Requisition  Size  (units) 

2 

2 

Average  Units  Demanded  per  Month 

.5 

2 

Average  Lead  Time  (Months) 

8 

8 

Coefficient  of  Variation  of  Lead  Time 

.25 

.25 

The  effects  of  varying  each  of  these  parameters  while  holding  the  other  three 
constant  were  investigated  for  base  cases.  We  found  that  the  SNB  approximation 
provides  percentage  points  very  similar  to  the  LPG  except  when  average  requisition 
sizes  exceed  k  units  or  when  demand  per  month  is  very  small.  Hence,  use  of  the  SNB 
approximation  appears  to  be  a  promising  way  to  significantly  reduce  the  computer 
requirements  associated  with  the  LPG  distribution. 

We  also  analyzed  the  skewness*  versus  kurtosis  relationships  for  the  sample  data 
sets.  These  relationships  indicate  that  the  continuous  Gamma  distribution  may  be  an 
excellent  approximation  to  the  Logarithmic-Poisson-Gamma. 
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Appendix  A 


Detailed  Item  Plots 


SAMPLE  LPG  CALCULATIONS 


A-  ] 


Data  Set  No. 

1 

_ 2 

_ 3 

_ 4 

_ 5 

Ax’e.  Req .  Si2e 

1.01 

2 

4 

8 

16 

Ave.  Demand/Mo. 

.5 

.5 

.5 

.5 

.5 

Mean  Lead  Time 

8.0 

8.0 

8.0 

8 

8 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 
Theta 

.02 

.72 

.90 

.96 

.99 

Lambda 

.49 

.25 

.12 

.06 

.03 

Alpha 

16 

16 

16 

16 

16 

Beta 

2 

2 

2 

2 

2 

LPG  Moments 

Mean 

4 

4 

4 

4 

4 

CV 

.56 

.97 

1.63 

2.63 

4.15 

Skewness 

.41 

1.46 

2.96 

5.08 

8.18 

Kurtosis 

3.23 

6.40 

16.45 

42.13 

103.79 

1 

Percentage  LPG  NB 

2 

rpG  Rs 

EPS — re 

lPg  RE 

5 

rpc — RB 

Points 


50 

4 

5 

3 

4 

1 

4 

0 

8 

0 

18 

60 

4 

5 

4 

5 

2 

5 

0 

10 

0 

22 

70 

5 

:  .  6 

5 

6 

4 

6 

2 

11 

0 

26 

80 

6 

7 

7 

7 

7 

10 

4 

13 

1 

29 

85 

6 

7 

8 

8 

9 

12 

7 

14 

3 

31 

90 

7 

8 

9 

9 

12 

14 

12 

22 

7 

33 

95 

8 

9 

12 

12 

17 

16 

22 

33 

22 

43 

97 

9 

10 

13 

13 

21 

22 

31 

38 

38 

68 

99 

10 

11 

17 

16 

31 

26 

52 

42 

82 

93 

Table 

II -2.  1 

Requisi 

tion 

Size 

Sensitivity 

for  LO 

Base 

Case 

DEMAND  IN  LEND  TINE  PARAMETERS 


EXACT  LPG  PROBABILITIES 
SCALED  NEGATIVE  BINONIAL 


A-6 


SAMPLE  LPG  CALCULATIONS 


Data  Set  No. 

_ 6 

_ 7 

_ 8 

_ 9 

10 

Ave.  Req.  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand/Mo 

.50 

.50 

.50 

.50 

.50 

Mean  Lead  Time 

4.00 

6.00 

8.00 

10.00 

12.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.25 

.25 

.25 

.25 

.25 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

4.00 

2.67 

2.00 

1.60 

1.33 

LPG  Moments 

Mean 

2.00 

3.00 

4.00 

5.00 

6.00 

Skewness 

2.16 

1.72 

1.46 

1.28 

1.14 

Kurtosis 

10.22 

7.67 

6.40 

5.65 

5.15 

6 

7 

8 

9 

10 

Percentage  LPG 

Points 

"HE  CFG 

Kb 

UFg  nb 

IPG 

TJB  LPG” 

NB 

.50  1 

2  2 

3 

3  4 

4 

5  5 

6 

.60  2 

2  3 

4 

4  5 

5 

6  6 

7 

.70  2 

;  4  4 

5 

5  6 

6 

7  8 

8 

.80  3 

5  5 

6 

7  7 

8 

9  9 

11 

.85  4 

5  6 

7 

8  8 

9 

9  11 

12 

.90  5 

6  7 

8 

9  9 

11 

12  12 

13 

.95  7 

8  10 

9 

12  12 

14 

14  15 

15 

.97  9 

9  11 

12 

13  13 

15 

15  17 

16 

.99  12 

9  15 

14 

17  16 

19 

19  22 

20 

Table  I I -3. 

Lead  Tine  Sensitivity  for 

LO  Base 

Case 

in  lend  Tine  pahaneteds 


EXACT  IPS  PROtAtlUTIES 
SCALED  NEGATIVE  DINOfllAL 


IN  LEAD  TINE 


—  EXACT  CPC  PAOIADILITIES 
-  SCALED  NEGATIVE  D1NOAIAL 


(MX  CPF  DIFFERENCE*  •.*•»  FOR  X 


--  EXACT  LPC  PROfAIIllTIES 
-  SCALED  NEGATIVE  IINONIAL 


A-12 


SAMPLE 

LPG  CALCULATIONS 

Data  Set  No. 

11 

12 

13 

14 

15 

Ave.  Req .  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand/Mo 

.10 

.50 

1.00 

1.50 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.05 

.25 

.50 

.75 

1.00 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2.00 

2.00 

2.00 

LPG  Moments 

Mean 

.80 

4.00 

8.00 

12.00 

16.00 

CV 

2.11 

.97 

.71 

.60 

.53 

Skewness 

3.52 

1.46 

.95 

.73 

.60 

Kurtosis 

21.73 

6.40 

4.53 

3.93 

3.64 

11 

12 

13 

14  15 

Percentage  LPG 

Points 

NB  LPG 

HjSgljfl 

.50 

0 

1  3 

4 

7  8 

11 

12  15 

15 

.60 

0 

2  4 

5 

8  9 

13 

13  17 

18 

.70 

1 

2  5 

6 

10  11 

15 

15  19 

20 

o 

00 

• 

1 

2  7 

7 

12  13 

17 

19  23 

23 

.85 

2 

2  8 

8 

14  14 

19 

20  26 

26 

.90 

3 

5  9 

9 

16  16 

22 

22  27 

28 

.95 

4 

5  12 

12 

19  19 

25 

26  32 

32 

.97 

5 

6  13 

13 

21  21 

28 

28  35 

35 

.99 

8 

8  17 

16 

26  25 

34 

33  41 

40 

Table  II -4.  Demand  Rate  Sensitivity  for  LO  Base  Case 


sIMNO  IN  LEND  TINE  NNNNNETENS 


EXACT  LAG  PtODADIUTIES 
SCALED  NEGATIVE  tINONIAL 


OCIMNO  IN  LEAD  TIRE  PAAAHETEW 


IH  LEAD  TINE  AARANETERS 


MX  CDF  DIFFERENCE 


A- 1 8 


SAMPLE  LPG  CALCULATION 


Data  Set 

No. 

16 

17 

18 

19 

20 

Ave.  Req. 

Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand /Mo. 

0.5 

0.5 

0.5 

0.5 

0.5 

Mean  Lead 

Time 

8 

8.0 

8.0 

8.0 

8.0 

CV  of  Lead  Time 

.01 

.25 

.50 

.75 

1.0 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

• 

.25 

.25 

.25  * 

.25 

.25 

Alpha 

1,000 

16.00 

4.00 

1.78 

1.00 

Beta 

1,250 

2.00 

.50 

.22 

.13 

LPG  Moments 

Mean 

4.00 

4.00 

4.00 

4.00 

4.00 

CV 

.94 

.97 

1.06 

1.20 

1.37 

Skewness 

1.61 

1.46 

1.21 

1.13 

1.29 

Kurtosis 

6.84 

6.40 

5.55 

5.15 

5.76 

Percentag 

16  17 

e  LPG  EB  TPC  TIB 

18 

efg — m 

19  20 

LPG  “KB  TBS  KB 

Points 

.50 

3 

4  3 

4 

3  4 

2 

3  2 

2 

.60 

4 

.  5  4 

5 

4  5 

3 

5  3 

4 

.70 

5 

'  6  5 

6 

5  6 

5 

6  5 

6 

,80 

7 

7  7 

7 

7  8 

7 

8  7 

8 

.85 

8 

8  8 

8 

8  9 

8 

9  9 

9 

.90 

9 

9  9 

9 

10  11 

10 

11  11 

12 

.95 

11 

12  12 

12 

12  13 

14 

14  15 

15 

.97 

13 

13  13 

13 

14  14 

16 

17  18 

19 

.99 

17 

15  17 

16 

19  19 

21 

21  25 

25 

Table  1 1-5.  Lead  Tine  Variability  Sensitivity  for  LO  Base  Case 
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SAMPLE  LPG  CALCULATIONS 


A-24 


Data  Set  No. 

21 

22 

23 

_ 24 

25 

Ave.  Req .  Size 

1.01 

2.00 

4.00 

8.00 

16.00 

Ave.  Demand /Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.02 

.72 

.90 

.96 

.99 

Lambda 

1.98 

1.00 

.  50^ 

.25 

.13 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2.00 

2.00 

2.00 

LPG  Moments 

Mean 

16.00 

16.00 

16.00 

16.00 

16.00 

CV 

.36 

.53 

.84 

1.33 

2.09 

Skewness 

.27 

.60 

1.35 

2.44 

4.02 

Kurtosis 

3.13 

3.64 

5.97 

12.32 

27.79 

Percentage  IK  ^NB  TO~TI5  IFifTB  IK  2~Tlg  IK"  N6 
roints  —  -  .  .  .  . 


.50 

IS 

16 

15 

15 

13 

14 

8 

13 

2 

22 

.60 

17 

18 

17 

18 

16 

19 

12 

18 

S 

26 

.70 

19 

19 

19 

20 

20 

23 

lfc 

26 

11 

31 

.80 

21 

21 

23 

23 

25 

27 

27 

35 

23 

37 

.85 

22 

22 

25 

26 

29 

31 

33 

39 

33 

56 

.90 

24 

24 

27 

28 

34 

35 

43 

47 

50 

75 

.95 

26 

27 

32 

32 

42 

42 

59 

61 

80 

93 

.97 

28 

28 

35 

34 

49 

46 

72 

67 

105 

101 

.99 

31 

32 

41 

40 

61 

57 

99 

88 

163 

148 

Tabic  II -6.  Requisition  Sire  Sensitivity  for  HI  Base  Case 
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MX  COT  DIFFERENCE- 


IN  UN)  TIME 


*  ^ 

SIS  1 
£f  I 
£ 

9m  £ 

IS  s 

?5  h 

t|J  o 


|  me  £  *5  3S 

C  I  »•  •«•»  - 


n3  s  E  .5  I 
:1*»  *  i»  Ss  l 

k  it  I  tils 


.•sssmsmssafssss::: 

^  i  i  i  i  i  it 

c5S8«S2£28B*S**S*55W 


ttS8S«8=£25RSR??5a35S? 


rwx  c tf  piFrE«HCC-  row  x 


SAMPLE  LPG  CALCULATIONS 


A-30 


Data  Set  No. 

26 

27 

28 

29 

30 

Ave.  Req  Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand /Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead  Time 

4.00 

6.00 

8.00 

10.00 

12.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

1.00 

1.00 

1.00 

1.00 

1  .00 

Alpha 

16.00 

16.00 

16.00 

16.00 

16.00 

Beta 

4.00 

2.67 

2.00 

1.60 

1.33 

LPG  Moments 

Mean 

8.00 

12.00 

16.00 

20.00 

24.00 

CV 

.71 

.60 

.53 

.49  .46 

Skewness 

.95 

.73 

.60 

.52 

.47 

Kurtosis 

4.53 

3.92 

3.64 

3.48 

3.38 

26 

27 

28 

29  30 

Percentage  LPG 

Points 

-"EE  TPS~ 

"TO 

LPG  51 

LPG 

"55  TFft 

55 

.50  7 

8  11 

12 

15  15 

19 

20  23 

23 

.60  8 

9  13 

13 

17  18 

21 

22  25 

26 

.70  10 

•11  15 

15 

19  20 

24- 

25  29 

29 

.80  12 

13  17 

19 

23  23 

28 

28  33 

33 

.85  14 

14  19 

20 

25  26 

30 

30  35 

36 

.90  16 

16  22 

22 

27  28 

33 

34  39 

39 

.95  19 

19  25 

26 

32  32 

38 

37  44 

44 

.97  21 

21  28 

28 

35  34 

41 

41  48 

48 

.99  26 

25  34 

33 

41  40 

48 

48  55 

55 

Table  II -7.  Lead  Time  Sensitivity  for  HI  Base  Case 


KBMftD  II*  UU  TI«e  <**»*«*£« 


PROBABILITIES 
ITIVE  BINOMIAL 


DENAND  IN  LEAD  Tint  P  MANE  TEW 


A-36 


SAMPLE 

LPG  CALCULATIONS 

Data  Set  No. 

31 

32 

33 

34 

35 

Ave.  Req  Size 

2.00 

2.00 

2 

.00 

2.00 

2.00 

Ave .  Demand /Mo. 

.10 

.50 

1 

.00 

1.50 

2.00 

Mean  Lead  Time 

8.00 

8.00 

8 

.00 

8.00 

8.00 

CV  of  Lead  Time 

.25 

.25 

.25 

.25 

.25 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

.05 

.25 

.50  . 

.  7S 

1.00 

Alpha 

16.00 

16.00 

16 

.00 

16.00 

16.00 

Beta 

2.00 

2.00 

2 

.00 

2.00 

2.00 

LPG  Moments 

Mean 

.80 

4.00 

8 

.00 

12.00 

16.00 

CV 

2.11 

.97 

.71 

.60 

.53 

Skewness 

3.52 

1.46 

• 

.95 

.73 

.60 

Kurtosos 

21.73 

6.40 

4 

.53 

3.92 

3.64 

31 

32 

33 

34  35 

Percentage  LPg 

Poi nts 

NB  TO 

NB 

CFC~ 

NB 

LPfi 

“NB  li*g 

“NB 

.50  0 

1  3 

4 

7 

8 

11 

12  15 

15 

.60  0 

2  4 

5 

8 

9 

13 

13  17 

18 

.70  1 

2  5 

6 

10 

11 

15 

15  19 

20 

.80  1 

2  7 

7 

12 

13 

17 

19  23 

23 

.85  2 

2  8 

8 

14 

14 

19 

20  25 

26 

.90  3 

5  1 

9  9 

16 

16 

22 

22  27 

28 

.95  4 

5  12 

12 

19 

19 

25 

26  32 

32 

.97  5 

6  13 

3 

21 

21 

28 

28  35 

34 

.99  8 

8  17 

16 

26 

25 

34 

33  41 

40 

Table  II  -8. 

Demand  Rate  Sensitivity  for 

HI  Base  Case 

1WWCMS 


DCIMN0  IM  LEAD  TIHC  MMKTW 
ICO  SIZE  DM 


Dcwmo  ii  lead  Tine  mkmctcm 


MX  CDF  DIFFERENCE* 


SAMPLE 

LPG  CALCULATIONS 

Date  Set  No. 

36 

37 

_ 38 

_ 39 

40 

Ave .  Req. 

Size 

2.00 

2.00 

2.00 

2.00 

2.00 

Ave.  Demand /Mo. 

2.00 

2.00 

2.00 

2.00 

2.00 

Mean  Lead 

Time 

8.00 

8.00 

8.00 

8.00 

8.00 

CV  of  Lead 

Time 

.01 

.25 

.50 

.75 

1.00 

LPG  Parameter 

Theta 

.72 

.72 

.72 

.72 

.72 

Lambda 

•• 

1.00 

1.00 

1.00 

1.00 

1.00 

Alpha 

10,000 

16.00 

4.00 

1.78 

1.00 

Beta 

1,250 

2.00 

.50 

.22 

.13 

LPG  Moments 

Mean 

16.00 

16.00 

16.00 

16.00 

16.00 

CV 

.47 

.53 

.69 

.88 

1.10 

Skewness 

.80 

.60 

.65 

1.03 

1.55 

Kurtosi s 

3.96 

3.64 

3.69 

4.86 

7.10 

36 

37 

38 

39  40 

Percentage 

Points 

We¬ 

51 

Lpg 

51 

IfrG  nB~ 

Lp6 

KB  IPS 

"55 

.SO 

is 

16 

15 

15 

14  14 

12 

13  10 

11 

.60 

17 

18 

17 

18 

17  18 

16 

16  14 

15 

.70 

19 

20 

19 

20 

20  20 

19 

20 

.80 

22 

22 

23 

23 

24  25 

26 

23  27 

27 

.85 

24 

25 

25 

26 

27  27 

30 

26  32 

33 

.90 

26 

27 

27 

28 

31  31 

35 

35  39 

40 

.95 

30 

29 

32 

32 

37  37 

44 

44  51 

51 

.97 

32 

32 

35 

34 

41  41 

50 

50  60 

61 

.99 

38 

36 

41 

40 

50  51 

64 

63  80 

80 

Table  17-9.  Lead  Tine  Variability  Sensitivity  for  HI  Base  Case 


MX  CDF  DIFFERENCE*  FOR  X 


IN  UNO  T«W 


Appendix  B 


LPG  Program  Source  Listings 


*vi 


B-l 


LIST 

10  REN  LP6.S— COMPUTE  EXACT  LP6  PRO!  USIN6  0/19/81  RECURSION 
20  REN 

30  BIN  T(200),T2(200) 

40  REN 

90  REN  TOO  -  T(X,K)  FOR  THE  CURRENT  X 

40  REN  T21K)  •  T<X-1,K>  FRON  PREVIOUS  X  CALCULATION 

70  FILES  LP 6.1 

80  SCRATCH  81 

90  REN  INPUT  PARAHETERS 

91  PRINT  "LPS.S DEBUS  *  PRINT  FIX)*  (1  OR  0)“ 

93  INPUT  18,  B7 

94  REN 

100  GOSUB  180 

110  REN  DO  LP6  RECURSION 

120  G0SUB840 

130REN  DO  SCALLEI  BINONIAL  CALCULATIONS 

140  80SUB  3180 

ISO  PRINT  "CONTINUE?  <Y  OR  N>" 

1SS  INPUT  At 

140  IF  At»"Y"  THEN  90 

170  STOP 

180  REN 

181  PRINT 

182  PRINT 

190  PRINT  "LP60.8 — EXACT  LPO  PROB  CALCULATIONS  USIN6  RECURSION" 

191  PRINT 

192  PRINT 
200  REN 

210  PRINT  "INPUT  AVE  REfi,  E(D),  E(LT),  C.OF.V  OF  LT  " 

220  INPUT  R0,B1 fE1 ,C 
230  URITE  N1,R0,B1,E1,C 

240  REN 

290  REN  SOLVE  FOR  THETA  -  T1 

240  60SUB  1710 
270  T1-8 

280  REN 

290  REN  ESTIHATE  PARAHETERS  FOR  LPO 

300  REN  S  •  STB  DEV,  B  »  VAR  TO  NEAN  RATIO 

310  8  •  C*E1 

320  B  «  8*8 /El 

330  Ut-  DI*E1 

340  REN 

390  B1-1/B 

340  At  »B1*E1 

370  R1—T1/1  <t-T1)  •LOOtt-TD  ) 

380  LI  »  U1/1E1*R1> 

390  REM 

400  PRINT  "AVE  RED  8I2E  *",81 
410  PRINT 

420  PRINT  "THETA  »",T1 , "LAMBDA  »",L1 

410  PRINT  "ALPHA-",  A1,"DETA  -\N1 

440  REN 


/ 


B-2 


430  REN 

400  N1»-T1/((1-T1  )*106<1-T1  )<»L1*A1/I1 ) 

470  PRINT 

410  PRINT  “HE AN  >%N1 
490  PRINT 

500REN  ESTIMATE  THE  FIRST  FOUR  NONENTS  Ml ,N2,N3,H4  OF  THELP6  RISTR1IUTI0N 
510  81  »  81*<1-TI ) 

520  C4  «  -L1*T1/<L0»<1.-TD) 

530  REN 
S40REH 

550  HI  «  A1*C4/D1 

SOOREN 

570  H2  *  A1*C4*(IHC4)/Dr2 

SSOREN 

590  H3  *  A1*C4*<(I1‘‘2)*(I+T1 )  +2*04*21/81*3 
OOOREN 

010  N4  •  11*3*(1+4*T1+T1*2)  ♦  B1*2*C4*(3.*A1  +  1.) 

020  N4  «  N4  ♦  0*lt*C4«C4 *A1  +C4*3*<3*A1+0. ) 

030  N4  •  A1*C4*N4/D1 *4 

040REN 

050  PRINT  "NONENTS" 

000  PRINT  “NT  N2  M3  H4" 

070  PRINT  HI ,H2,M3,H4 

080REN 

090REN  COMPUTE  STANIARDIZEB  NONENTS 

700  PRINT 

710  S  -  S8R1H2) 

720  PRINT  “CQEF  OF  VAR  **,  S/N1 

730  PRINT  "N3/S*3«  ",  N3/S*3 

740  PRINT  "N4/S*4»  ",  N4/S*4 

750REN 

700REN  OUTPUT  LP6  PARAMETERS  TO  FILE 

770  REN 

780  NRITE  II ,M1 ,S/Nt ,H3/S*3,N4/S*4 

790  WRITE  II ,T1 ,L1,A1 ,11 

800  RETURN 

81  OREM  •♦•••♦•♦•♦•♦••♦•♦•••♦•••♦M****************************** 

820REN  LP6  RECURSION  CALCULATION 

830 

840  REN 

850  REN  REDEFINE  C  FOR  PROIAIILITY  CALCULATIONS 

BOO  REN 

870  C>  -L1/L08O  .-T1 ) 

880  Cl  >  C/<L1  +  ID 

890  REN 

900  REN  SET  LINIT  ON  *  ■  240 

910  REN 

920  L2*200 

930  REH 

940  C2*0 

950  REN 

900  REN  SET  CONSTANTS  FOR  USE  IN  RECURSION 


B-3 


970  REN 
990  S7»t 

990  H1«<B1/<L1*B1>rA1 

1000  H1-S7*H1 

1010  REN  S7  >  SCALE  FACTOR 

1020  N2-<C/<L1«B1)) 

1030  REN 

1040  PRINTa>C,C1,H1,H2M,C;C1;Hl;N2 
1050  REN 

1000  REN  EVALUATE  H(X«0) 

1070  REN 

1010  T(0)  •  0 

1090  T2(0)  *  HI 

1100  St  >  HI 

1110  C2  »  C2  +  S1/S7 

1120  REN 

1130  PRINT 

1131  URITEN1 r0,S1 ,C2 

1140  IF  B7<*  0  THEN  1200 

1141  PRINT 

1150  PRINT  *  X  H(X)  F(X)“ 

1151  PRINT 

1100  PRINT  o;si;c2 

1170  REN 
1110  REN 

1190  REN  EVALUATE  H(X)  FOR  X  >  0 
1200  REN 

1210  FOR  X  >  1  TO  L2 

1220  81-0 

1230  T<0)*0 

1240  REN 

1250  T<X)-(T1*H2*<A1*X-1 )/X)*T2(X-1) 

1200  REN 

1270  IF  X  <  2  GOTO  1380 

1210  REN 
1290  REN 
1300  REN 

1310  FOR  K-1  TO  X-1 

1320  T(K)-(T1/X)  •  <H2*<A1«K-1)*T2(K-t>  ♦<X-1)*T2(K>  ) 

1330  81-SHT1K) 

1340  IF  B8  <■  0  THEN  1300 

1350  PRINT”X,KvT<X(K),81->*,X;K;T<K>,’S1 

1300  NEXT  K 

1370  REN  PICK  UP  T(X,X>  TERN  IN  SUN 

1380  81-81  ♦  T(X) 

1390  REN 
1400  REN 
1410  REN 


PRINT  TOTALS  FOR  H<X> 


B- 


1420  C2  •  C2  ♦  S1/S7 

1430  PRINT  X;St;C2 

1440  MITE  N1  ,X,S1  fC2 

14S0  REN 

1460  REN  IF  CUN  PROS  EXCEEBS  .09,  STOP 

1470  REN 

1480  IF  C2  >  .99  60T0  1670 

1490  REN 

1500  REN  PRINT  T(K)  TERNS  FOR  DEBUGS INS 

1510  REN 

1520  IF  D8  <■  0  THEN  1560 

1530  FOR  K»0  TO  X 

1540  PRINT  *X,K,T(X,K)  *>*,X;X;T(K) 

1550  NEXT  N 
1500  REN 
1570  REN 

1580  REN  RECORB  TOO  VALUES  FOR  USE  IN  NEXT  PASS 

1590  REN 

1600  FOR  K-0  TO  X 

1610  T2<K)«T(K> 

1620  NEXT  K 

1630  REN 

1640  REN - END  OF  X  LOOP 

1650  NEXT  X 

1660  REN 

1670  WRITE  11,-99,-99,-99  •  END  OF  LPO* 

1680  RETURN 

1690  REN»****>»<m»***»m*«**««**'**********»«**REND  OF  LPO  RECURSION 
1700  REN 
1710  REN 
1720  REN 

1730  REN - AVE  REQ  SIZE  EBUATION 

1740  BEF  FNR(fl)  >  -8/(  (1-0)*L06(1-0)  ) 

1750  REN 

1760  REN - 

1770  REN  BINARY  SEARCH  ROUTINE 

1780  REN - - - 

1790  REN 

1800  REN  SXVE  FOR  THE  VALUE  OF  Q  UNICH  8IVES  AND 
1810  REN  AVE.  RE8UISI7I0N  SIZE  OF  RO 
1820  REN 

1830  REN  SET  UP  END  POINTS  FOR  SEARCH 

1840  REN 

1850  X-  0 

1860  Q9-.999 

1870  R9>FNR(Q9) 

1880  81>.001 
1890  R1«FNR(81 ) 

1900  SOTO  2080 
1910  REN 


CHECK  IF  RO  <«  R 


1920  REN 
1930  REN 
1940  IF  RO  >  R  THEN  2000 
1930  REN 

1940  REH  RE-SET  TOP  OF  INTERVAL 

1970  R9«R 

1980  Q9>0 

1990  GOTO  2060 

2000  REN 

2010  REN  RE-SET  BOTTON  OF  INTERVAL 

2020  R1«R 

2030  01 *0 

2040  SOTO  2080 

2030  REN 

2040  REH  HALVE  THE  INTERVAL,  AND  RE-EVALUATE  THE  FUNCTION 

2070  REN 
2080  *«*♦? 

2090  Q«(0t*Q9>/2 

2100  R-FNR(Q) 

2120  REN 

2130  REN  IF  R  IS  VITKIN  .001  OF  RO,  THEN  RETURN 

2140  REN 

2130  IF  ARS(R-RO)  <  .001  THEN  2110 

2140  IF  K  >  15  THEN  2180 

2170  80T0  1940 

2180  RETURN 

2190  REN 

2200REN  LPO.S— LP  PROSAIIUTV  CALCULATIONS 

2210REN 

2220REN  USELPG.B  AS  THE  OUTPUT  FILE 

2230REH 

2240  REN 

2230  SCRATCH  81 

2240  PRINT  -THETA,  LAHBIA,  ALPHA,  SETA,  DEBUG?" 

2270  IHPUT  T1 , L  t ,A1 ,11 ,18 

2280  PRINT  "EXACT  LP6  PROS* 

2290  605UB  2420 

2300  PRINT  "SCALEI  NE8.  BIN” 

2310  BOSUI  3230 

2320  PRINT 

2330  PRINT  ”CONTINUEW  OR  N>” 

2340  INPUT  At 

2330  IF  At>”Y"  GOTO  2240 

2340  STOP 


B-6 


2370  REN  IP6.S--HASTER  LPO  ROUTINE  LIBRARY 
2380  REH  LP61.S 

23 tO  PRINT  "THIS  PROGRAM  CONFUTES  EXACT  PRO!  FOR  THE  LFG  BIST* 
2400  BIN  T ( 1 00, 1 00) 

2410  INPUT  T1,L1,A1,B1 

2420  PRINT  “THETA  >“;T1tMLANIDA  »“;L1 “ALPHA  »“$A1 ,MBETA  »";B1 
2430  URITE  11,1,2,3,4 

2440  URITE  81 ,Tt ,L1 ,A1 ,B1 

2450  C  ■  -  LI  /LOGO  -T1) 

24A0  HI  ■  -T1/<(l-T1 )*L06(1-T1  D*L1*At/B1 

2479  Cl  >C/(L1  ♦  ID 

2480  L2  «  INT(I00*N1 ) 

2490  C2  -0 

2SOO  FOR  X  -0  TO  12 
2S10  Y(X,0>-0 
2520  FI  -1 

2330  IF  X<2  THEN  2570 
2540  FOR  K  ■  1  TO  X-1 
2550  FI  •  F1*K 
2500  NEXT  K 
2570  Y  <X,X)  •  1 ./FI 
2580  FI  >  1 

2590  IF  XO  0  THEN  2430 

2400  SI  -<B1/(L1  +I1))*A1 

2410  60  TO  2740 

2420  REN 

2430  SI  >0 

2440  FOR  K  -  1  TO  X 

2450  IF  X  -  1  THEN  2470 

2440  Y(X,K)*Y(X-1 ,K-1 )/<X-1 )  4  Y<X-1,K> 

2470  FI  >  I 
2480  FOR  J  -  0  TO  K-1 
2490  FI  •  F t * < At  4  J14C1 
2700  NEXT  J 

2710  81  -  St  4  Y(X,K)*F1 
2720  NEXT  K 

2730  SI  ■  S1*(81/ai  4  B1D'A1*T1*X/X 

2740  C2  ■  C2  4  81 

2750  IF  08  <»0  THEN  2780 

2740  PRINT  X,S1,C2 

2770  WRITE  81,X,S1,C2 

2780  IF  C2  >  .99  THEN  2800 

2790  NEXT  X 

2800  URITE  01,-99,-99,-99 
2810  RETURN 
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2820  REN 

2830  REN  MMI*t****M*»l*MMM«**Mt*t«Mt**Mt*»*UMt««t*»M 

2840  REN  LP62.S 

2130  PRINT  "SCALES  POISSON" 

2800  INPUT  T1,L1,A1,B1 

2870  PRINT  "THETA  »"JT1 , “LAMIA  «"JL1, "ALPHA  «";A1,"BETA  «";B1 
2880  C  ■  -L1/L06<Wn 
2890  C3»  1./(I.-T1) 

2900  Kt  -  T1*C 

2910  C4  ■  (B1/(K1+81 ))*A1 

2920  U  «  KI/(Kf481) 

2930  C2  *  0 

2940  FOR  N  «  0  TO  1000 

2950  PI  •  1 

2900  IF  NOO  THEN  3000 
297#  PI  ■  C4 

2980  SOTO  .3040 

2990  REN 

3000  FOR  I  «  0  TO  N-1 

3010  PI  •  P1*( (A1+N-I-1 )/<H-I )*U) 

3020  NEXT  I 

3030  PI  »  P1*C4 

3040  C2  -  C2  ♦  PI 

3050  IF  B8<>0  THEN  3080 

3000  PRINT  N;C3*N;Pi;C2 

3070  UR1TE  81 fC3*H,P1 ,C2 

3080  IF  C2>.99  THEN  3100 

3090  NEXT  N 

3100  RETURN 

3110  RETURN 

3120  REN  •**••*•••«*•«***•*«****•*«***•******»*•*•••*»*••*»•»*»**** 

3130  REN 
3140  REN  LP63.S 

3130  PRINT  "THIS  PR08RAN  CONFUTES  APP  PROS  FOR  LPS  DIST  US INS  SCALES  SIN" 

3100  INPUT  Tl yLt t A1 fit 

3101  PRINT 

3170  PRINT  "THETA  «";T1 ,"LAHBDA  -";11,"ALPNA  -";A1,"BETA  »";BI 
3180  PRINT 

3190  PRINT  "NEC  SINONIAL  PROSAIILITIES" 

3200  PRINT 

3210  PRINT  "X  P  F" 

3220  PRINT 

3230  C  •  1/(1. -Tl) 

3240  XI  •  -T1*LI/L0C(1-T1) 

3250  C4  >  (B1/(K1 +B1 ) ) *A1 
3200  U  •  K1/4K1+B1 ) 

3270  C2  •  0 

3280  FOR  N  -  0  TO  1000 
329#  PI  ■  1 

3300  IF  NOO  THEN  3370 
3110  PI  ■  C4 
3320  A  •  -1 

3330  B  -  INTtC/2  *0.5) 

3140  82  ■  P1/(C/2*l ) 

3130  SOTO  3440 
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33*0  REN 

3370  FOR  I  ■  0  TO  N-1 

3380  PI  •  P I  *  < ( A1 )/<N-I)*U) 

3390  NEXT  I 

3400  PI  ■  P1*C4 

34T0  4  «  INT<(2*H-1 >*C/2+0.5) 

3420  I  ■  INT ( <2*N+1 >*C/2+0 .5 ) 

3430  S2  •  PI/C 
3440  L3  ■  C2 
3430  FOR  X  >  4+1  TO  I 
3440  F3  ■  C2  ♦  S2*<X-4) 

3470  F2  -  F3  -  L3 
3480  L3  »  F3 

3490  IF  87  <»  0  THEN  3520 

3500  PRINT  X,F2,F3 

3510  URITEN1 ,X,F2,F3 

3520  NEXT  X 

3530  C2  -  C2  +-  PI 

3540  IF  C2  >  .9?  THEN  3540 

3550  NEXT  N 

3540  WRITER 1,-99, -99, -99, "ENB  OF  NB" 

3570  RETURN 
3580  REN 

3590  REN  •***«******•***•*«**•***»*»«**»*«»* 

3400  REN  LP64.S 

3410  REN  COHPUTE  REQUISITION  SIZE  R  VS  THET4  T4BLE 

3420  REN 

3430  FOR  T1  «  .01  TO  .99  STEP  .01 
3440  8  «  —  < 1  — T 1 )*L06(1-T1) 

3450  R  «  T1/8 
3440  PRINT  T1,R 
3470  NEXT  T1 
3410  STOP 
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3490  REN 

3700REH  DEFINITIONS  IN  NAHNIAS'  P 
3710  REN  DEFINITIO 
3720  REN 

3730  REN 

LAND 

LI 

3740  REN 

ALPH 

Al 

3750  REN 

BET 

Bt 

3740  REN 

NEAN 

HI 

3770  REN 

CNST 

Cl 

3780  REN 

UN 

12 

3790  REN 

CUN 

C2 

3800  REN 

FACT 

FI 

3810  REN 

SUN 

SI 

3820  REN 

KN 

K1 

3830  REN 

PROD 

PI 

3840  REN 

CN 

C3 

3850  REN 

CONST 

C4 

3840  REN 

LAST 

L3 

3870  REN 

FX 

F2 

3880  REN 

FFX 

F3 

3890  REN 

3900  END 

3L 

82 

